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Abstract 

Metal-insulator transitions driven by disorder (A) and/or by electron correlations (U) are investigated within the 
Anderson-Hubbard model with local binary-alloy disorder using a simple but consistent mean-field approach. The 
A-U phase diagram is derived and discussed for T = and finite temperatures. 
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If spontaneous symmetry breaking is excluded, a 
system of electrons in a non-degenerate half-filled va- 
lence band may undergo a transition from a normal 
Fermi liquid to an insulator either due to Coulomb 
interaction or due to disorder. Metal-insulator transi- 
tions (MIT) in the presence of strong electron correla- 
tions and disorder are not well understood - even on 
the mean-field level. For the purely correlation-induced 
(Mott) MIT, the dynamical mean-field theory (DMFT) 
has uncovered a rather complex phase diagram [1]. 
The MIT to an (Anderson) insulator in case of non- 
interacting electrons and strong diagonal binary-alloy 
disorder can be described by the coherent-potential ap- 
proximation (CPA) [2]. While spatial correlations are 
neglected in both cases, the residual mean-field physics 
at low temperatures T is non-trivial and relevant for 
three-dimensional transition-metal oxides, for exam- 
ple. The combined problem can be studied within the 
half-filled (n = 1) Anderson-Hubbard model (AHM): 
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Here the n.n. hopping is set to t = 1, U is the on-site 
interaction, fj, — U/2 the chemical potential, and e; = 
±A/2 with equal probabilities x = 1/2 a random on- 
site energy at site i. A measures the disorder strength. 
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Fig. 1. U-A phase diagram for T — (insets: U-T phase 
diagram for different A). Energy scale: t — 1. 

"DMFT-I-CPA" [3] can be regarded as the optimum 
mean- field approach to this model. This, however, must 
be supplemented by stochastic [3] or renormalization- 
group techniques [4] or by further approximations, e.g. 
weak-coupling perturbation theory [5]. While parts of 
the phase diagram in the U-A-T space are known [3,5], 
a comprehensive study is still missing. 

Here we employ the self-energy-functional approach 
(SFA) [6] and the Us = 2-site dynamical-impurity 
approximation (2S-DIA) which nicely reproduces the 
phase diagram for A — [6]. In case of disorder, 
a proper generalization of the formalism has to be 
applied [7]. Within the generalized framework, the 
2S-DIA can be regarded as a strongly simplified but 
consistent DMFT-CPA approach. In the limit ris — » cx) 
one recovers the DMFT for A = 0, the CPA for [/ = 
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and the DMFT-CPA for (7, A / 0. 

Operationally, the Green's function G' of a single- 
impurity Anderson model H' with two sites and im- 
purity on-site energies e = ±A is obtained by exact 
diagonalization and averaged, T' — {G'), to get the 
configuration- independent self-energy S' = Gq^ — 
r' where Go is the free {U, A = 0) Green's function. 
S' — S{t') depends on the one-particle parameters of 
H' and is used as a trial self-energy in a general varia- 
tional principle, (5f2[S] = 0, which gives the exact av- 
eraged grand potential of the AHM at the physical S. 
On the subspace given by S{t'), the functional can be 
evaluated rigorously (see [7]). We consider the para- 
magnetic phase of the AHM on a three-dimensional 
s.c. lattice consisting of 10"' sites. Phase boundaries 
are obtained from the resulting f2 as a function of A, 
U and T. The averaged interacting local density of 
states (DOS) of the AHM can be calculated via p{ij) = 
-lmrii{uj + iri)/n and T = (G,7^ - S')'^ where Go is 
the free ([/, A = 0) lattice Green's function. 

Three different phases are identified at T = (see 
Fig. 1): a paramagnetic metallic phase (PM), a Mott 
insulator (MI), and an Anderson insulator (AI). For 
any disorder strength A, we find the AI at weak U 
(and A > Ac(Lf)) to be well separated from the MI at 
strong U by the PM in between. For A = the criti- 
cal interaction for the Mott MIT is found to be f/c = 
13.9 ^ 1.161^ (with W ^ 12 the free band width) 
while Ac = 5.4 = QAQW for the MIT at (7 = 0. This 
agrees well with full DMFT and CPA estimates, re- 
spectively [1,2,6]. For f/ci < U < Uc2, a coexistence 
of the stable PM phase with the metastable MI phase 
is observed (f/d = 12.4). This scenario for the Mott 
MIT is well known for A = [1] and is shown here 
to survive for any finite disorder strength with a A 
dependent coexistence region f/ci(A) < U < Uc2{^) 
and t/c(A) = [/c2(A). A discontinuous Mott MIT with 
f7ci(A) < Uc{A) < [/c2(A) is found for finite tempera- 
tures < T < rc(A). For T > Tc(A) there is a smooth 
crossover only. For A — > cx), the critical interactions ap- 
proach a linear dependence, [/ci,2(A) ^ A -|- consti,2 
while T'c(A) -^ Tc saturates. 

The topology of the phase diagram can be under- 
stood by looking at the DOS, see Fig. 2. Characteris- 
tic for the MI at A = and U > Uc is the insulating 
gap between the lower and upper Hubbard band (LHB, 
UHB). For finite A the gap decreases due to a broad- 
ening and, eventually, a splitting of each of the Hub- 
bard bands (Fig. 2, A = 6). The closure of the gap is 
preempted by the occurrence of a quasi-particle peak 
(QP) at u> — which marks the transition to the PM 
(see A = 14). Apart from the QP, the spectrum can be 
understood as being composed of two Hubbard bands 
at oj ~ ±A/2 ± [7/2 for each of the two atomic config- 
urations e = ±A/2. This explains the strong spectral- 
weight transfer when increasing A = 18 to A = 20: 
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Fig. 2. Average density of states for (7 — 18 and different A 
and Lorentzian broadening with rj — 0.25. 

As the e = -I-A/2-DOS (e = -A/2-D0S) becomes al- 
most completely unoccupied (occupied), the weight of 
the UHB (LHB) must disappear. Finally, a further in- 
crease of A induces a splitting into an upper and lower 
alloy band (UAB, LAB) and a MIT to the AI. 

Recently, Byczuk et al [4] have shown that 
DMFT+CPA predicts the AHM to exhibit a Mott 
MIT also for fillings n 7^ 1 if a:: = n. Similar to the 
presently considered case n = 1 = 2x, a sharp QP 
at w = (even for strong disorder) as well as a co- 
existence of the PM and the MI is found close to the 
MIT. We like to point out that the phase diagram for 
n = X =^ 1 [4] can be understood by an analysis of the 
DOS completely analogous to the n = 1 = 2a:: case dis- 
cussed above - although its topology is quite different. 

Concluding, we have proposed a mean-field scenario 
for the MIT in the AHM at half-fiUing n = 1 = 2a: on 
the basis of a simplified DMFT-I-CPA approach. The 
phase diagram can be understood by a quasi-atomic in- 
terpretation of the DOS in combination with the Mott 
MIT scenario of the pure system. This should be con- 
trasted with full DMFT-I-CPA calculations in the fu- 
ture which may also clarify the importance of disorder 
scattering due to a finite self-energy ImS{iJ = 0) which 
has been neglected here. 
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